Ground-state energy density of a dilute Bose gas in the canonical transformation 
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A ground state energy density of an interacting dilute Bose gas system is studied in the canonical 
transformation scheme. It is shown that the transformation scheme enables us to calculate a higher 
order correction of order no? in the particle depletion and ground state energy density of a dilute 
Bose gas system, which corresponds to the density fluctuation contribution from the excited states. 
Considering two-body interaction only, the coefficient of no? term is shown to be 2(ir — 8/3) for the 
particle depletion, and 16(n — 8/3) for the ground state energy density. 
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I. INTRODUCTION 

Investigation of an interacting dilute Bose gas has a long history beginning with the seminal works of Bogoliubov, 
Beliaev, Lee, Huang and Yang, Hugenholtz and Pines. This topic has been amply discussed in standard text booksi 2 -^ 
Despite the successful demonstration of Bose-Einstein condensation (BEC) in inhomogeneous system of dilute alkali 
gases during the last decade, the dilute interacting Bose fluids in homogeneous system is one of the fundamental 
interests still in statistical physics. 

The ground state energy density of the dilute system is expressed in the following general form^ 

f —EE^'^^F, (i) 

i=0 3=0 

where n = N/V is the particle density. Since unconfined BEC occurs for a two-body repulsive interaction, the s- 
wave scattering length a is positive, [i/2] is the Gauss number which takes the integer part of i/2, and are the 
coefficients in the expansion. Through elaborate calculations, it is known thal£ 
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Comparing Eq. Q with Eq. ©, we find the coefficients as C 00 = 1, C w = 128/150?, and C 2 i = 8(4tt - 3y/3)/3. 
Nevertheless, the coefficient of C21 was calculated early in 1950'sj^i>£ the constant k = C20 and other higher order 
terms did not receive much attention. 

Hugenholtz and Pines mentioned that the expansion is not a simple power series, rather it involves the logarithm of 
expansion parameter na? ? Although Fetter and Walecka gave a rough estimate on the fluctuation contribution from 
the condensate in the order of no?, they stated that the coefficient C20 has never been determined^ In order to get 
a deeper understanding of the ground state over the mean field results, it is necessary to calculate corrections from 
the quantum fluctuations. Earlier derivation of rigorous upper and lower bounds had been given by Dyson4 Also an 
improved calculation on the lower bounds has been carried by Lieb and Yngvason.- 

Some practical calculations of the coefficient C20 were carried out very recently by two groups. Our group obtained 
C20 from the fluctuation contribution by using the functional Schrodinger picture approach of two-body interaction^ 
The Schrodinger picture for many particle systems or relativistic quantum field theory is an infinite dimensional 
. extension of quantum mechanics and leads into a functional formulation. 11 i 12 i 13 However, the result has not been 
confirmed using a separate independent method. Another contribution by three-body interaction has been calculated 
by Braaten et al&^ In their calculation the core potential is necessary to find C20, and they calculated C20 as a 
function of an effective potential parameter r s £ Therefore, the two results for C20 are independent. 

In this paper, we consider the contribution from the condensate fluctuation, which is independent on the details 
of the potential and two-body interaction only. We show that the result obtained by the Schrodinger picture can be 
obtained by a more standard canonical transformation method. In the process we obtain the fluctuation correction 
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to the particle depletion, and Cio in Eq. ([I]) from second order ground state expansion. Also, we will consider the 
fluctuation contribution from the condensate, and then calculate the coefficients C20 and C30. On the other hand C21 
which come from three-body interaction will not be touched. 

This paper is organized as follows: In Section II, the interacting Hamiltonian of the dilute Bose system is expanded 
up to the terms of our interest. In Section III, the canonical transformations are applied twice to diagonalize the 
Hamiltonian. In Section IV, the particle depletion appearing in the diagonalized Hamiltonian is obtained. In Section 
V, the fluctuation correction to the ground state energy density is obtained as a series of no 3 . In Section VI, we 
summarise and discuss our results. 



II. MANY-BODY DILUTE BOSE GAS SYSTEM 



We start with the standard model Hamiltonian for an interacting Bose gas given by2. 
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Here, the constant matrix element g should be determined by requiring that H produces the two-body scattering 
properties in vacuum. Since the system is dilute enough, we assume that the s-wave scattering length a is much less 
than the inter-particle spacing n 1 / 3 . As a result of the na 3 <C 1, this satisfies to be a small expansion parameter. 
Under this assumption the zero momentum operators ao, a J become nearly classical due to occurrence of a condensate. 

In two-body interaction the interacting part of the Hamiltonian in Eq. ([3]) is expanded to keep the terms up to 
order of V 2 , N , and //Vq as±&±£ 
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The last term which is generally neglected in textbook calculations originates from the interactions of particles out 
of and into the condensate, thus representing the condensation fluctuation. In Ref. [5] the three-body interaction 
term (ip*^) 3 contributes to the C2o- However, in our calculation one momentum among the kj., k2, k3, k.4 of two-body 
interaction is zero and will contribute to the C2o- Then, it makes our contribution different from that of Ref. [5]. 
In order to handle the contribution from the last term, we introduce a new variable jk 
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Replacing the operators ao, a by VVo> and applying the number operator relation Vo = N — X]fc^o( a I afc + a -fc a -<0/2 
to Eq. we rewrite the model Hamiltonian as 

H = \n 2 Vg 
2 y 
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where e9 = h 2 k 2 /2m. The ground state properties of the above Hamiltonian is to be calculated using a canonical 
transformation method. 
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III. CANONICAL TRANSFORMATION METHOD 



To diagonalize the above Hamiltonian, we apply the Bogoliubov transformation of the new the operators a k and 



a\ as follows ;i 



a 
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Then, the Hamiltonian in Eq. ([6]) becomes 



(bl+A k b_ k ). (7) 
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is chosen to make the off-diagonal term vanish. That is, 

2(e° k +ng)A k + ng(l + A 2 k )=0, (9) 



and we obtain A k {= A- k ) as 



where 



Ak = ^lM±M, (10) 
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^ (e° + ng) 2 - (rig) 



E k = J(& + na)*-(na)». (11) 



The positive sign was selected in Eq. (JTUJ) to give positive excited energy values. Note that (1 + A k ) / (1 — A k ) — e^/E k . 
Substituting A k into the Hamiltonian in Eq. ([8]), we obtain the following compact form of the Hamiltonian 

H = U 2 Vg + J2r^A2 [{*° k +ng)Al + ngA k ] 

1 fc#0 k 

+ Y,E k b\b k + Y,G k (b k +bl k ), (12) 

fc#0 k^O 

where 

We can regard the system as a collection of quantum mechanical harmonic oscillators exposed to additional forces 
given by linear terms. The linear terms b k and b\. can be made to eliminate by a linear transformation of the form 

b k = c k + a k , 

b\ = 4 + a_ fc) (14) 

where the new operators c k and ct satisfy the bosonic commutator relations. Note that at = a_ k and commutes 
with Cfc. 
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Substituting the new operators in Eq. JT4J) into Eq. (fl2|) . we can rewrite the Hanhltonian in Eq. (|12|) as 

H = \n 2 Vg + T-^-H M + n 9) A l + n 9 A k] 

1 k^O 1 A k 

+ ^2 (E k c\c k + E k a\a k + 2G k c 



k^Q 



+ Y^{E k a- k + G k )c k + Y^{Eka k + G-k)4' (15) 

k^O k=iO 

Here, we choose a k = —G- k /E k to make the linear terms vanish. Then, the Hamiltonian in Eq. (|15[) is now 
diagonalized in terms of new operators as follows 



H = \n 2 Vg + J2 Y^tf K e ° + n 9) A l + n 9 A k] 



fc#0 

G k G_ k 



E k 

fe#0 k^O 

= Hg + ^EkUk- (16) 

Using Eqs. (flQ]) . (fTT|) . (fl"3]) . and the chosen ak, the above Hamiltonian ff s produces the general form of the ground 
state energy as 

E g = \n 2 Vg + i]T [E k - (e° + ng)] + ^ £ ^< l7fc7-fc |0) 

= £ + £i + £ 2 . (17) 
From the relation of the interaction strength, g,— 
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we obtain the ground state energy density of the first two terms as 
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Note that J2 k is replaced with V J d 3 k/(2ir) 3 during the calculation. The first two terms, Eq and E\ in the ground 
state energy are exactly same as those obtained through conventional methods in the literature 

IV. THE PARTICLE DEPLETION 

Accurate determination of the fractional depletion of the zero momentum state for an interacting gas is increasingly 
needed with the advance of the experimental technique on the BEC states. The particle distribution of the dilute 
system is obtained from Eq. Q and (fT4l as 
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The dominant part of the unknown quantity (■y k ~f^ k ) is obtained by iteration method. Applying the definition of 
7^ in Eq. ([5]), it is written as 

<0| 7fe 7- fe |0) = <0| Yl aU q a q al k+ql a q ,\Q). (21) 
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It is known that the dominant contribution arises when a particle interacts with itself and it belongs to the term 
where q' = k + q. Then, we separate the dominant contribution into two parts q = q' and q^q' again. 
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+ (k dependent minor terms). 

The term for q ^ q' vanishes using the argument of random phase approximation. 
Within this approximation we obtain the dominant part of ("fkj-k) as 
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The new variable p = ^/na 3 /ir (0 < p <C 1) accounts for the diluteness of the system and will be used for the 
expansion. 

When we back {jk'Y-k) m Eq. (f2"3"| into Eq. (|2H)) by iteration method, we obtain a new nk as 
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Putting the new rik into Eq. (|23p . we obtain {ykj-k) as 
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We see the dominant part of ("fkl-k) is fc-independent. Therefore, E 2 in Eq. (|17p can be written as 
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Finally substituting Eq. (|25|) into Eq. (|20D . we obtain the particle depletion from the zero momentum condensate 
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We obtained the next order correction term of the particle depletion beyond the textbook result as 2(ir— 8/3)na 3 !^^ 
If we repeat the above procedure, we can obtain the higher order terms too within the given approximation scheme. 
In the following section, we discuss the next order correction to the ground state energy density which originates from 
the fluctuation contribution of the condensate. 



6 



V. THE HIGHER ORDER CORRECTION TO THE GROUND STATE ENERGY DENSITY 

This approximation of ("fk"/-k} also gives us a simple result for the condensation fluctuation correction E2 in Eq. 
(f2"6")) to the ground state energy density, but the integral carries an ultraviolet divergence. The same divergence in 
Ei was handled through the expression of the scattering length a to order of g 2 as given in Eq. (TT5)) . However, for 
the calculation of E% such a simple cancellation scheme is not possible, because it is now necessary to expand a up to 
g 3 . Therefore, we resort to an alternative cutoff procedure called minimal subtraction renormalization scheme. The 
justification for this scheme to the current problem is given in detail in Ref. [5]. 

In this scheme the ultraviolet divergence in Eq. (|26[) can be regularized by imposing a cutoff procedure in the k 
space. With the cutoff |k| < A where A is an arbitrary separation between short distance effect, the divergence can 
be handled by introducing a counter term in the Hamiltonian. That is, linear, quadratic, and other power ultraviolet 
divergences are removed as part of the regularization scheme by subtracting the appropriate power of k from the 
momentum space integrated. In this case 1/fc 2 is subtracted from the integral. The justification for this procedure 
is that the terms that are subtracted are dominated by short distances and can be canceled by counter terms in the 
Hamiltonian. Also, we note that the result is independent of the potential used, since the fluctuation contribution 
from the condensate considered here is a generic feature for any potential. 3 

Following this prescription, we obtain Ei in a dimensionless form 
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Substituting Eqs. (TT8"|) . (f!H)) and (129j) into Eq. (l2"51) . this approximation gives the result for the condensation 
fluctuation correction, E%, 

g = 16. (tt - |) p 2 + 32. (n " |) (y " *) P 3 + Otf). (30) 

Thus, the present calculation gives C20 = 16(tt — 8/3), and C30 = (32 / y/ir) (ir — 8/3)(10/3 — ir). The C20 is exactly the 
same as one from our previous functional Schrodinger picture method.— The higher order contribution to the ground 
state energy in Eq. (|30[) comes from the interactions of particles out of and into the condensate. This term represents 
the fluctuation contribution from the excited states and disappears if the number of condensate particles is conserved 
strictly. 



VI. SUMMARY AND DISCUSSIONS 



In this paper we revived the ground state energy density calculating of an i nteracting dilute Bose gas. Hugenholtz 
and Pines have given a general form of the series expansion in terms of V na 3 . Till now only the first few coefficients 
are known. In view of the fundamental interest of homogeneous dilute interacting Bose gas, it is desirable to know 
about the coefficients of higher order expansion. 

There are at least two possible sources to give contributions for the unknown coefficients of na 3 terms. One is the 
fluctuation from condensation and produces a universal result. The other is from the three-body interaction and not 
universal because it depends on the radial potential. Here we used the standard canonical transformation method 
of two-body interaction. We obtained the next order terms in the particle depletion and the ground state energy 
density coming from the condensate fluctuation. They are 2(ir — 8/3) and 16(7r — 8/3) respectively. Furthermore, the 
coefficient of (na 3 ) 3 / 2 term is (32/ v / ^)( 7 r - 8/3)(10/3 — tt) in this approximation. 
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